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ON ; Abstract 

We show that the GBW saturation model provides an exact solution to the one- 



dimensional linear transport equation. We also show that it is motivated by the BK 



Qh' equation considered in the saturated regime when the diffusion and the splitting 

Q^' term in the diffusive approximation are balanced by the nonlinear term. 

> 

Perturbative Quantum Chromodynamics (pQCD) at high energies can be formulated in 
. coordinate space in the dipole picture []]. If we in particular focus on Deep Inelastic 



1 Introduction 



Scattering the scattering process can be described in this picture as interaction of virtual 
photon which has just enough energy to dissociate into a 'color dipole' with the hadronic 
target carrying most of the total energy. The interaction process is described here by 
the dipole-nucleus scattering amplitude. The formal formulation and application of the 



dipole picture leads to high energy pQCD evolution equations for the dipole amplitude 
and in particular to a BFKL equation for this quantity [2] which together with a certain 
factorisation theorem allows to calculate observables in the dipole picture. However, leav- 
ing out the formal approach one can also model the dipole amplitude. The Golec-Biernat 
Wiisthoff saturation model developed roughly 10 years ego [3] is a quite well tested model 
for the dipole amplitude which includes saturation effects. It was motivated by require- 
ments that at the high energy limit of QCD the total cross section for hadronic processes 
should obey unitarity requirements. At present there are much more sophisticated ap- 
proaches to introduce these requirements in a description of scatterings at high energies 
[U El El E E] • However, one can still ask the question if there is any dynamics behind the 
GBW model or to put it differently is there any equation to which formula proposed by 
Golec-Biernat and Wiisthoff is a solution? And what is the role of the initial conditions? 
In this article we want to answer these questions. 

The Letter is organized as follows: in the next section we show that the GBW model 
provides exact solution to the one-dimensional transport equation. We show this for the 
unintegrated gluon density and for the dipole amplitude in the momentum space. In the 
third section we relate the derived transport equation to the BK equation which allows 
us to conclude that the GBW model is well based in pQCD. 



2 GBW model and a transport equation 



2.1 Momentum space analysis 

The GBW amplitude following form GBW cross section and related to it by a(x, r) = 
2 J d 2 bN(x, r, b) reads (here we are interested in the original formulation without evolution 
in the hard scale [9]): 



N(x,r,b) = 9(b -b) 



1 — exp 



2 



4R% 
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where b is the impact parameter of the collision defined as distance between center of the 
proton with radius bo and center of a dipole scattering on it, r is a transversal size of the 

dipole, x is the Bjorken variable, Ro(x) = ^ ( ^ j is the so called saturation radius and 
its inverse defines saturation scale, Q s (x) = l/Ro(x) and xo, A are free parameters. This 
amplitude saturates for large dipoles r ^> 2R Q and exhibits geometrical scaling which has 
been confirmed by data [TO] . 

2.1.1 Transport equation for unintegrated gluon density 

The dipole amplitude ([TJ can be related to the unintegrated gluon density which convo- 
luted with the kx dependent off-shell matrix elements allows to calculate observables in 
the high energy limit of QCD. This relation is the following (after assumption that the 
dipole is much smaller than the target) [TTJ [T2"] : 



t( 72 n N c ,4v72 f d2r , -i ,N(x,r,b) 
f(x, k ,b) = -k V k / — exp -zk • r 



(2) 



where r and k are two-dimensional vectors in transversal plane of the collision and r = |r| 
k = |k| 

Performing this transformation we obtain the known result 



f(x, k\b) = .r^-eibo - b)R 2 (x)k 4 ex V [-R 2 (x)k 2 ] (3) 

Now motivated by the fact that this formula exhibits a maximum both as a function of x 
for fixed k 2 and as a function of k 2 for fixed x, we differentiate f(x, k 2 , b) with respect to 
x and f(x, k 2 , b)/k 2 with respect to k 2 . We obtain: 

W (^, l) = VWfip« (4) 

xQq 



a k 2 U TJTxi ^ 



f( x,k 2 ,b) = f(x,k 2 ,b)(l-R 2 (x)k 2 ) 
k 2 WQl 
Dividing eqn. (jl]) by ([5]) and rearranging the terms and defining T[x, k 2 ,b) = f(x, k 2 ,b)/k 2 , 
Y = lnxo/x, L — In k 2 jQ\ we obtain: 



d y T{Y, L, b) + \d L T{Y, L, b) = (6) 



which is the first order linear wave equation also known as the transport equation. As it 
is linear it cannot generate saturation dynamically but it can propagate well the initial 
condition leading to a successful phenomenology [3] . It describes the change (wave) in the 
particle distribution flowing into and out of the phase space volume with velocity A. This 
wave propagates in one direction. The quantity T{x, k 2 , b) gains here the interpretation 
of a number density of gluons with momentum fraction x with the transversal momentum 
k 2 at distance b from the center of the proton. The general solution of ([6]) can be found 
by the method of characteristics and is given by: 

F(Y,L,b)=F (L-\Y,b) (7) 

One can go back from to Q using following initial condition at x = Xq\ 

F(x = x ,k 2 ,b) = -^e(b -b)k 2 exp(-k 2 ) (8) 

This initial condition has saturation built in, since the gluon density vanishes for small 
k 2 . One can also try a different initial condition at large k 2 : 

F(x = x ,k 2 ,b)~^ (9) 

which gives solution without the saturation effect: 

Knowing the properties of the linear first order partial differential equation we see that the 
property of saturation of GBW was a consequence of the wave solution which relates x and 
k 2 supplemented by initial conditions with saturation built in. We also see that the critical 
line of the GBW saturation model visualizing, the dependence of the saturation scale on 

A/2 

x, Q s (x) = Qo (—J is in fact from the mathematical point of view the characteristics 
of the transport equation. 



2.1.2 Transport equation for the dipole amplitude in momentum space 

Similar investigations can be repeated for the momentum space representation of the 
dipole amplitude N(x, r, b) which we denote by <p(x, k 2 ,b). 

, 2n f d2r i -i ,N(x,r,b) , . 

(f)(x,k 2 ,b)= I — exp(-zk-r) ^' 2 ' (11) 

A nonlinear pQCD evolution equation like the Balitksy-Kovchegov (BK) equation written 
for <p (in large target approximation) takes quite simple form and can be related directly 
to the statistical formulation of the high energy limit of QCD (see [H] and references 
therein). Applying this transformation to ([T|) we obtain: 

^{x ) k 2 1 b)= l -6{b -b)T 



(12) 



where T(0, z) = J 2 °° ye~ 2 . To obtain the partial differential equation we proceed similarly 
as before and upon differentiating (Tl2|) with respect to k 2 and also with respect to x and 
introducing Y = lnxo/x, L = Ink 2 jk\ we obtain: 

d Y (f)(Y, L, b) + \d L <f>(Y, L,b) = (13) 



which is, as before, the transport equation which supplemented with the initial condition, 



{x = x ,k 2 ,b) = -e(b -b)T 



k 2 



(14) 



gives back (Ti"2j) . 



2.2 Relation to pQCD 

It is tempting to investigate the relation between ([6]) and the high energy pQCD evolution 
equations like [51 EJ [8]. Let us focus here in particular on the form of the BK equation in 
large cylindrical target approximation for the dipole amplitude in momentum space for 
which the nonlinear term is just a simple local quadratic expression. The BK equation 
for the dipole amplitude in the momentum space reads: 

d Y <j>(Y, k 2 , b) = a X (-^—) <t>(Y, k 2 , b) - a^iY, k 2 , b) (15) 

V °log k 2 J 

where a = and x{l) = 2^(1) — ~ rfO- ~ l) is the characteristic function of 
the BFKL kernel which allows for emission of dipoles and therefore drives the rise of the 
amplitude. The role of the nonlinear term is roughly to allow for multiple scatterings 
of dipoles which contributes with negative sign and slows down the rise of the ampli- 
tude. This equation provides unitarization of the dipole amplitude [T5] for fixed impact 
parameter and admits traveling wave solution in the diffusion approximation [T6] . 



2.2.1 Analytic approach 

The analytic solution of ffT5|) within the diffusion approximation relying on expanding the 
kernel of ( Tl5l) up to second order and mapping it to the Fisher-Kolmogorov equation has 
been obtained by Munier and Peschanski [T7]. It reads: 



k 2 \ ( k 2 ^ 7C_1 



In 2 



A- 2 



2a X "(lc)Y \Q 2 (Y)J\ 
(16) 

where 7 C = 0.373 and Q 2 (Y) is emergent saturation scale given by: 



By inspection we see that (flBT) does not obey the transport equation. The problem is 
caused by the diffusion term. However, we can consider the asymptotic regime called 



"front interior" [TTJ Q2|, region where transverse momenta k is close to the saturation 
scale Q S (Y) and rapidity Y is large and where the condition In 2 ^ q^jy) ) /^x"(lc)Y 
is satisfied. In this regime ([To]) simplifies to: 



«y, *>, *) = - »)^- ) ( JgJ 7 " loglf/cKOOl (18) 

Proceeding as in the previous sections we obtain the following wave equation: 

dY<P(Y, L, b) + X BK d L <P(Y, L, b) = (19) 

where Xbk = d log Q 2 S (Y) /dY . In the limit where Xbk does not depends on energy [Jj5] 
we obtain: 

Xbk = -ax'ilc) (20) 

2.2.2 Numerical results 

Now let us investigate the exact numerical solution of the BK equation in the diffusion 
approximation to see how its terms arrange to reduce to the transport equation in the 
saturation regime. Following [16J we represent the BFKL kernel of the BK equation as a 
power series around ■y c : 



xhd L )4>(Y,L,b) 



X(7c) + (~d L - 7c)x'(Tc) + j x {-d L - 7c ) Y'(7c) + - 



<f>(X,L,b) 

(21) 

Taking terms up to the second order and rearranging them we obtain the following equa- 
tion: 

d Y <P(Y,L,b) _ _ D<p(Y,L,b)-<j ) \Y,L,b) 

~ a , (v t u\ = a Xohc) ~ a (22) 

d L (j){Y, L, b) d L (f){Y, L, b) 

where D = Xoilc) ~ IcXoilc) + \(r d L ~ lc) 2 Xo(lc)- 

The right hand side of this equation using the analogy to the transport equation should 
define velocity of the wave. Solving numerically ([22]) with the initial condition 0(0, L, b) = 
9(bo — b)e~ L2 we obtain results for wave velocities which are shown on Fig. (JT^.). 
Defining ratio: 

D<j>{Y,L,b)-4> 2 {Y,L,b) 

X ' ( lc )d L c!>(Y,L,b) { ] 

and plotting it we observe (Fig. ([lb)) that in the saturation region (small L, large Y) 
the second term in ( l22|) is proportional to the first one and approaches 2«Xo(7c) which 
together with the first term gives the velocity of a wave traveling towards higher values 
of log A; 2 . Using this fact we can write in that in asymptotic limit Y — >oc: 

" W(YTM =- q *o(7 c ) = Ab* (24) 



which is the same as (|T9l) . The numerical value Xbk = 0.92 is clearly different as compared 
to GBW approach where velocity A = 0.27 (four flavor fit) is a free parameter which is to 
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Figure 1: Velocity of wave plotted as a function of Ink 2 for four different values of rapidity: 
Y=5, 10, 20, 30 (a). Ratio R defined in the text for four different values of rapidity: Y=5, 10, 
20, 30. This plot shows that at saturation region the BK equation reduces to transport equation. 

be determined by data. However, in the case of the BK equation A is a number only in the 
asymptotic region. When contributions beyond asymptotic ones are taken into account 
like diffusion, proper kinematics in emission of dipoles, renormalization group effects [20] 
it becomes a quite complicated function of energy. 

3 Conclusions 

In this note we have shown that the GBW saturation model is the exact solution of a 
one-dimensional linear transport equation of the form (JHD- We conclude that since ($5$) 
is a linear equation the saturation property has to be provided in the initial condition. 
We found that for the GBW model this equation is universal for the unintegrated gluon 
density f(x, k 2 , b) and the dipole amplitude in momentum space <p(x, k 2 , b) but the details 
of the shape of the wave depends on the initial condition which is different for each of 
them. We also studied the relation of the transport equation to the BK equation in the 
diffusion approximation. We have shown that in the region of phase space where diffusion 
and splitting processes are of the same order as the nonlinear term, the GBW model is 
consistent with the BK equation. 
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